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ABSTRACT 

We present new gas flow models for the Milky Way inside the solar circle. We use SPH 
simulations in gravitational potentials determined from the NIR luminosity distribu- 
tion of the bulge and disk, assuming constant NIR mass-to-light ratio, with an outer 
halo added in some cases. The luminosity models are based on the COBE/DIRBE 
maps and on clump giant star counts in several bulge fields, and include a spiral arm 
model for the disk. 

Gas flows in models which include massive spiral arms clearly match the observed 
^^CO diagram better than if the potential does not include spiral structure. 

Furthermore, models in which the luminous mass distribution and the gravitational 
potential of the Milky Way have four spiral arms are better fits to the observed {l,v) 
diagram than two-armed models. 

Besides single pattern speed models we investigate models with separate pattern 
speeds for the bar and spiral arms. The most important difference is that in the latter 
case the gas spiral arms go through the bar corotation region, keeping the gas aligned 
with the arms there. In the (Z, v) plot this results in characteristic regions which appear 
to be nearly void of gas. In single pattern speed models these regions are filled with 
gas because the spiral arms dissolve in the bar corotation region. 

Comparing with the ^^CO data we find evidence for separate pattern speeds in 
the Milky Way. From a series of models the preferred range for the bar pattern speed 
is rip = 60 ± 5Gyr~^, corresponding to corotation at 3.4 ± 0.3 kpc. The spiral pattern 
speed is less well constrained, but our preferred value is Vlsp « 20Gyr^^. A further 
series of gas models is computed for different bar angles, using separately determined 
luminosity models and gravitational potentials in each case. We find acceptable gas 
models for 20° < ^bar < 25°. The model with ((/^bar = 20°, % = 60Gyr"i, Vl^p = 
20Gyr~^) gives an excellent fit to the spiral arm ridges in the observed (/, v) plot. 

Key words: Galaxy: structure, Galaxy: kinematics and dynamics. Galaxy: centre, 
Galaxies: spiral. Interstellar medium: kinematics and dynamics. Hydrodynamics. 



1 INTRODUCTION 

Observations of cold gas in the Milky Way [MW] have con- 
tributed substantially to our understanding of MW struc- 
ture. No other tracer is observed in as large a part of the MW 
as are gas clouds . Lon gitude-velocity [Iv] d iagrams (Hart- 
mann & Burton (1997), Dame et al. (2001)) show the dis- 



tribution of gas velocities as a function of galactic longitude 
I, integrated over some range in latitude b. By observing the 
MW in different spectral lines, this gas can be traced at sub- 
stantially different densities. The largest absolute velocity as 
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a function of I defines the terminal velocity curve (TVC). In 
an axisymmetric galaxy, the gas at these velocities is found 
at the "tangent point" where the line of sight is tangential to 
a circle around the Galactic centre. From this the rotation 
curve can be determined. However, due to the bar and spiral 
perturbations in the MW potential, the gas has substantial 
non-circular velocities, which are most evident in the central 
10° - 20°, due to the bar, but also as "bumps" in the TVC 
where spiral arm tangents perturb the gas flow by ~ 10- 
20kms~^. At sub- TVC velocities, crowding in both position 
and in velocity produces ridge-like structures in the {l,v) 
diagram. The Galactic spiral arms are visible as straight or 
curved such ridges. 

A number of attempts have been made to model these 



© 0000 RAS 



2 Nicolai Bissantz, Peter Englmaier, and Ortwin Gerhard 



observations. One group is formed by analytic models of spi- 
ral structure. The first exhaustive analytical for mulati on of 
a spiral arm theory was developed b y Lin fc Shu ( 1964 ) and 
applied to the MW by Lin et al. (1969). They proposed 
a two-armed model with pitch angle —6° and a pattern 
speed for the spiral structure Qsp ~ 13.5km/s/kpc. Ama- 
ral & Lepine (1997) fitted the rotation curve of the MW to 



t ern s peeds in their A^'-body simulations. Rautiainen & Salo 
( 1999| ) analysed two-dimensional A'^-body simulations, some 



an analytic mass model and found a self-consistent solution 
with a combi ned t wo- and four-armed spiral structure. In 
Lepine et al. ( 2001 ) they extended this model to allow for a 
phase difference between the two-armed and the four-armed 
spiral pattern. 

The second group, numerical simulations of the Galac- 
tic gas flow, also have a long tradition. A recent example 
is the smoothed particles hydrodynamics (SPH) models of 
Fux (1999), who evolved a gas disk inside a self-consistent 
AT-body model scaled to the COBE/DIRBE K-band map of 
the MW and the radial velocity dispersion of M giants in 
Baade's window. The resulting gas flow was transient, but 
at specific times resembled closely a number of observed 
arms and clumps in the bar region. Weiner & Sellwood 
(1999) compared predictions from fluid dynamical simula- 



tions for analytic mass densities with the observed outer 
velocity contours of the HI {I, v) diagram to constra in th e 
pattern speed and bar angle. Englmaier & Gerhard (199!:) 



[hereafter. Paper I] computed gas flows in the gravitational 
potential of the NIR l uminosity distribution of Binney, Ger- 
hard & Spergel ( 1997 ), assuming constant NIR mass-to- light 
ratio (M/L). Their best SPH gas flow models reproduced 
quantitatively a number of observed gas flow features, in- 
cluding the positions of the flve main spiral arm tangents at 
|/| < 60° and much of the terminal velocity curve. 

An important feature of all these gas flow models is 
the pattern speed of the non-axisymmetric component. En- 
glmaier & Gerhard found a best pattern speed for the 
bar ilp « 60km/s/kpc. Weiner & Sellwood derived a bar 
pattern speed « 42km/s/kpc, whereas Fux d etermined 
« 50km/s/kpc from his models. Dehnen (2000) used res- 



onant features in the Hipparcos stellar velocity distribu- 
tion to argue that the Sun is located just outside the OLR 
of the exciting quadrupole perturbation, giving a pattern 
speed of ~ 51km/s/kpc for solar constants -Ro = 8kpc 
and Vo = 220 kms_^_. For their spiral stru cture model, 
Amaral & Lepine ( 1997 ) and Lepine et al . (2001 ) found 
fisp « 20 Gyr"' -35 Gyr"\ Fernandez et al. (|200l| ) obtained 



a somewhat higher D,sp ~ 30km/s/kpc from Hippa rcos data 
for OB stars and Cepheids. Debattista et al. (2002) used the 
Tremaine- Weinberg method on a sample of intermediate age 
to 8 Gyr old OH/IR-stars in the inner Galactic disk. They 
found a pattern speed of 59±5± 10 (systematic) km/s/kpc, 
which may be driven by the bar in the center of the MW. 

Thus the bar and spiral arms in the MW may not ro- 
tate with the same pattern speed. For a fast bar, a single 
ilp would imply that the spiral arms are entirely outside 
their corotation radius. Observations of external g alaxie s 
(see, e.g., the Hubble Atlas of Galaxies, Sandage (1961)) 



suggest that galaxies exist with dust lanes on the inner 
(concave) edges of their spiral arms. For a trailing spiral 
pattern, these arms would be inside their corotation ra- 
dius. A lower pattern speed for the spiral structure than 
for the bar woul d rem ove this discrepancy. Indeed, Sell- 
wood & Sparke ( 1988 ) showed evidence for multiple pat- 



of them with a massless, dissipative gas component added. 
They confirmed the possibility of multiple pattern speeds in 
self-consistent N-body models of barred galaxies, and found 
a number of possible configurations. These included mod- 
els with corotating bar and spirals, as well as models with 
different pattern speeds. Some of the models in the latter 
group sho w evi dence for a non-linear mode-coupling (Tag- 
ger et al. (1987)) between bar and spiral pattern, but others 
show no such evidence. In some of their models there ex- 
ist separate inner spirals corotating with the bar, and outer 
spirals which rotate with their own, lower pattern speed. It 
is therefore tempting to analyse gas fiow models of the MW 
with multiple pattern speeds. 

What is the morphology of the MW spiral arms? Most 
authors infer four spiral arms from tracers which directly 
or indirectly measure the gas density, such as molecu- 
lar clouds, HII regions, pulsa rs, an d the galactic magnetic 
field (Georgelin & Georgelin (|l97^), Sanders et al . (1985|), 



Caswell & Hayne s ( [19871 ), Grab elsky & et al. ( |1988|), Ta ylor 
& Cordes ( |l993| ), Vallee ( |l995| ); however. Bash ( |1981D in- 
fers a two-armed pattern from the same Hll-data as used by 
Georgelin & Georgelin). The problem is that all spiral arm 
parameters other than the tangent point directions (e.g.. 
Paper I) require distance information. It is also not clear 
whether all spiral arms seen in the MW gas are present 
in the old disk. Ortiz & Lepine ( 1993| ) constructed a four- 
armed model which reproduces their star counts in the near 
infrared. Drimmel (200C) preferred a two-armed structure 
from COBE/DIRBE K-band data, but a four-armed struc- 
ture for the du st dist ribution seen in the 240^m data. Drim- 
mel & Spergel ( 2001 ) project a luminosity model through the 



240/im dust model, to compare with the NIR J and K band 
COBE/DIRBE data. Their best model for the stellar dis- 
tribution is four-armed, but dominated by two spiral arms. 
Drimmel & Spergel conclude that, if there are four arms in 
the K-band luminosity distribution, the Sag-Car arm is of 
reduced strength (by a factor of 2.5). 

In this paper we investigate the dynamical effects of 
the Galactic bar and spiral arms on the gas flow in the 
Milky Way. We investigate the possibility of different pat- 
tern speeds for bar and spiral arms, and the consequences 
this would have on the observed {I, v) diagrams. Our mass 
models for the inner Galaxy are based on t he NI R luminos- 
ity density models of Bissantz & Gerhard (2002) [hereafter. 
Paper II] which include spiral structure. We use SPH sim- 
ulations to determine the gas flow in the MW, extending 
the work presented in Paper I, where eightfold symmetric 
mass models were used. This paper is organised as follows. 
In Section ^ we describe the luminosity models, methods, 
and observational data used in this work. In Section ^ we 
describe our best gas model for the observed ^^CO {l,v) 
diagram. Then we compare models with different pattern 
speeds (Section ^), bar angles and spiral arm morphology 
(Section ^) to constrain these parameters, and finally give 
our conclusions in Section M. 



2 GAS DYNAMICAL MODEL 
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2.1 Mass model of the Milky Way 

Here we first describe the adopted model for tlie distribu- 
tion of luminous mass in tlie MW. Tfie model is based on the 



dust-corrected near-infrared maps of Spergel et al. (1995) 



which they obtained from COBE/DIRBE data using a three- 
dimensional dust model derived from 240pim observations. 
From their L-band map we obtained a non-parametric lumi- 
nosity distribution using the procedure described in Bissantz 
& Gerhard (2002). Because the non-parametric model only 



covers the central 5kpc part of the MW, we used a para- 
metric best-fit model of the same L-band data to extend the 
model to larger radii. To convert the luminosity model to a 
model for the luminous mass density we assumed constant 
L-band mass-to-light ratio. 

The model building procedure is described in detail 
in Paper IL Here we only summarise the most important 
points. Bissantz & Gerhard estimated the luminosity model 
iteratively from the L-band data by maximising a penalised 
likelihood function. The penalty terms encourage eightfold- 
symmetry with respect to the three main planes of the 
bar, as well as smoothness and a prescribed spiral structure 
model in the disk. The spiral structure term is based on an 
approximate model for the MW spiral arms from Ortiz & 
Lepine ( [l993| ). The best models resulting from this approach 
reproduced the dust-corrected COBE/DIRBE maps with an 
RMS accuracy of 0.07 mag. 

To define the eightfold symmetry penalty term we had 
to specify the position of the Sun in the MW. Bissantz & 
Gerhard set the distance of the Sun to the galactic centre 
to 7?o = 8 kpc, and the vertical distance from the galactic 
plane to zo = 14 pc. These values will be used throughout 
the present paper as well. The third parameter is the bar 
angle. Paper H compared photometric models for bar angles 
10° < ^bai < 44°, and concluded that the best models had 
bar angles 20° < ^bar < 25°. 

Paper H also showed that models which include spiral 
structure are better than models without spiral structure. 
In particular, because some spiral arm tangent points are 
evident in the L-band data, such models give a better de- 
scription of the nearby disk. Also, the inclusion of spiral arms 
makes a model appear broader on the sky. For given data, 
this allows the bulge/bar to be more elongated in models 
with spiral structure. The larger bulge elongation in these 
models makes it possible to reproduce the asymmetries seen 
in the apparent magnitude distributions of clu mp gian t stars 
in several bulge fields (Stanek et al. (1994), (1997)), for 



15° ^ ifibar ^ 30°. The shape of the bulge/bar in the model 



with <^bar ~ 20° is about 10 : 3 
approximately 3.5 kpc. 



4:3 — 4, and its length is 



2.2 Gravitational potential 

The gravitational force field and potential generated by the 
distribution of luminous mass were calculated using the mul- 
tipole expansion method described in Paper I, modified to 
include phase terms in order to allow for the spiral arm com- 
ponents. The multipole expansion method was used because 
(i) it aUows an independent treatment of disk, bar, and spi- 
ral arm components, (ii) it works for an arbitrary density 
distribution, and (iii) it does not require detailed bound- 
ary conditions as do other methods, since the integrations 



Figure 1. Rotation curve of the standard mass model for bar 
angle 20° . The velocities have been scaled with the factor ^ de- 
termined in Section || fitting the observed terminal velocities by 
the SPH model in 10° < |«| < 50°. 



Figure 2. Quadrupole and octopole terms of the standard poten- 
tial, separated into the part generated by the triaxial bar/bulge 
in the distribution of luminous mass (full line: m = 2, dotted line: 
m = 4), and the part caused by the spiral arms (short dashed: 
m = 2, long dashed: m = 4). All multipoles are normalised by the 
value of the monopole term at the respective Galactic radius. 



can easily be extended to infinite radii. The expansion was 
computed to maximum spherical harmonic orders l,nax = 8 
and rrimax ~ 8. Odd orders were neglected, assuming point- 
symmetry with respect to the centre. Naively, one could in- 
terpret the m = 2 component as due to the bar mode and 
all higher m modes as due to spiral arm modes. However, 
some of the m — 2 component beyond the bar's corotation 
radius is due to an m = 2 spiral ar m mo de in the luminosity 
model (see also Amaral & Lepine (1997)). 



Since the spatial resolution of the mass models obtained 
in Paper H is limited by that of the dust-corrected maps of 
Spergel et al., the central cusp in the MW's density distribu- 
tion and potential is incorrectly represented in our multipole 
expansion model. We attempted to correct for this by mod- 
ifying the multipole coefficient functions in the centre as in 
Paper I, replacing the central mass distribution by a power 
law p"^ **. 

Optionally we add an analytical halo potential 
0Haio = i'K'fln(r'-f a"), 

where Vi^t is the circular velocity at infinity and a is the core 
radius. In this paper we use Vi^t = 220km/s and determine 
the core radius such that the model in the potential with 
halo best fits the TVC near \l\ = 90°. 

Figure |l] shows the rotation curve derived from the 
m=0-component of the reference luminosity density model 
of Paper II (for bar angle 20°, with the parametric exten- 
sion). For brevity, this mass density is called "standard mass 
model" hereafter, and the associated gravitational poten- 
tial "standard potential". Its rotation curve no longer shows 
a strong bump in the inner Galaxy's rotation curve as for 
the models of Paper I (Fig. 4 there). This is because the 
new luminosity model reproduces light in the MW's spiral 
arms significantly better, as compared to the density max- 
ima ~ 3 kpc down the minor axis of the bar in the earlier 
models. We note that the spiral arms in the luminosity mod- 
els used here arbitrarily start at a galactocentric radius of 
3.5 kpc. This generates a slight distortion of the rotation 
curve around this radius. 

Figure ^ shows the quadrupole and octopole terms of 
the potential, separated into the parts generated by the bar 
and the parts generated by the spiral arms. All multipoles 
are normalised by the value of the monopole term at the 
respective Galactic radius. Note that in this model the bar 
plays a major role for the non-axisymmetric forces, in par- 
ticular the quadrupole moment, even beyond where it ends 
in the mass density. 

The multipole representation of the potential is used 
in the orbital analysis and in the two-dimensional SPH hy- 
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Figure 3. Resonance diagram for the standard mass model, as- 
suming the same velocity scaling as in Fig. ^ 



drodynamics code, as described in Paper I. The bar and 
spiral components are both given a constant pattern speed. 
In most models described below, the bar and spiral arm 
patterns rotate with different pattern speeds, implying that 
the mass distribution and potential undergo periodic oscil- 
lations; in some models these pattern speeds are equal and 
the mass distribution is constant in the rotating frame. 

Figure ^ shows a resonance diagram for the standard 
mass model with the same scaling as in Figure |l[ (The scal- 
ing constant is different by < 1.5% for models 40 and 60, 
which have significantly different spiral pattern speeds from 
our standard model.) We will see below that the preferred 
bar pattern speed is Sip — 60Gyr"\ Observe from Fig. ^that 
in this case the corotation radius of the bar nearly coincides 
with the inner ultra-harmonic resonance of the spiral struc- 
ture, if this has pattern speed flap = 40Gyr"^, and with its 
inner Lindblad resonance, if Qsp ~ 20Gyr"^. Models with 
these values for the spiral arm pattern speed in conjunction 
with Qp — 60 Gyr"^ are particularly interesting to analyse 
(Tagger et al. (1987)), and are discussed below. 



2.3 Hydrodynamical method 

We use the two-dimensional smoothed particles hydrody- 
namics [SPH] code described by Englmaier & Gerhard 
(1997). The code solves Euler's equation for an isothermal 
gas with effective sound speed Ca: 



V)v 



(1) 



This approach is based on a result of Cowie (1980). He 



showed that an isothermal single fluid description crudely 
approximates the dynamics of the ISM. However, here the 
isothermal sound speed is not the thermal sound speed, but 
an effective sound speed representing the RMS random ve- 
locity of the cloud ensemble. 

The SPH method has the advantage of allowing for 
a spatially adaptive resolution length. This is achieved by 
adjusting the smoothing length /i of a particle everywhere 
such that the number of particles that overlap a given parti- 
cle is approximately constant. Fluid quantities are approx- 
imated by averaging over neighbouring particles. Further- 
more, the SPH scheme includes an artificial viscosity to allow 
for shocks in the simulated gas flow. For fu rther discussion 
of the method see Engl maier fc Gerhard (1997), Paper I, 
and Steinmetz & MiiUer ( |l993| ). 

Our SPH models contain 50000 — 60000 particles, except 
when indicated otherwise. The initial surface density of the 
models is taken to be constant inside 8 kpc galactocentric 
radius. 



2.4 Comparison with observational data 

Our main tools to compare our gas flow models with obser- 
vations are the terminal velocity curve (TVC) and the {I, v) 
diagram, which shows the radial velocities of gas clouds as 
function of galactic longitude. Throughout this paper ob- 
served velocities are given with respect to the local stan- 



dard of rest (LSR). For the distance of the Sun from the 
galactic centre we assume -Ro = 8 kpc. The LSR circular 
velocity is assumed to be Vg = 220km/s, consisten t wit h 
Ro = 8 kpc (Feast fc W hitelock ( |l997| ), Reid et al. ( [l999| ). 
Backer & Sramek (1999)). Note that a 10% change in vq is 



not critical, amounting to radial velocity variations of only 
~ lOkms^^ in the central I = 45° (cf. Paper I). 

{l,v) diagrams for our gas flow models are constructed 
as follows. The LSR observer is specified by the galactocen- 
tric radius Ro, the LSR circular velocity Vo, and by the angle 
95bar relative to the bar. We first project all particle veloc- 
ities {vx,Vy) onto the line-of-sight vector Cp from the LSR 
observer to the particle, subtracting the component of Vo in 
the direction of the particle: 



Vr = Cp ■ (vx, Vy) — Vo sin(/). 



(2) 



We then construct a two-dimensional binned histogram of 
the particle distribution in the / — w-plane, with bin size 
« 0.23 X 0.7km/s. Finally, we convert this histogram into 
a greyscale plot, using a lower surface density cutoff Cw ~ 
0.5 — 1% to enhance the contrast. Cw varies between different 
models, and is selected so as to optimize the visibility of the 
spiral arm ridges and the terminal velocity envelope. 

We compare a model with observations in a two step 
process. In the first step, we focus on the terminal veloci- 
ties, comparing the model TVC with an observed TVC com- 
posed from the following data: HI velocities from Burton & 
Liszt (1993); HI velocit ies fr om Fich et al. (1989), based on 
data from Westerhout (1957); unpublished 140-ft single dish 



HI velocities, kindly provided by D r. B. Burton; northern 
^^CO-velocities from Clemens ( 1985| ), including error bars; 
and southernjf CO-velocities from Alvarez et al. (1990). The 
Clemens ( 1985 ) data were corrected for internal dispersion 
(Paper I). The observed velocities, corrected by the respec- 
tive authors or Paper I to the pre-Hipparcos LSR frame, 
in which the Sun was assumed to move with approximately 
uq — — lOkms"^ inwards and Vq = 15kms~^ in the for- 
ward direction of Galactic rotation, are here transformed to 
the Hipparcos LSR frame {uq = — 10kms~\ Vq = 5kms^^, 
Dehnen & Binney (199S)), by subtracting 10 sin Z km s^^. 

A free parameter of our models is the mass-to-light ra- 
tio. From the comparison of the model TVC with the ob- 
served TVC we determine the best-fit scaling factor ^ for 
the model velocities. This parameter ^, which determines 
the mass scale of the model, varies by ~ 5 — 10% between all 
our models. We only use the TVC data for 10° < \l\ < 50° 
to determine ^, because near the centre the resolution of 
our models is insufficient, and because for \l\ > 50° the halo 
contributes significantly (see below). In the fitting we take 
special care of the location of "bumps" in the TVC, because 
these indicate spiral arm tangents. This procedure assumes 
that the NIR disk and bulge are responsible for all of the 
observed velocities in the central parts of the MW. 

In most models the gravitational potential is time de- 
pendent, because the pattern speeds of bar and spirals are 
different. In this case we select a "best" snapshot, corre- 
sponding to a specific phase and evolutionary age. The value 
of ^ generally depends slightly on both the model and the 
evolutionary age of the snapshot. 

In the second step we compare the model {l,v) dia- 
gra m wit h the observed ^^CO {l,v) diagram of Dame et 
al. (^00 1|) (Figure 0), using the scaling ^ from the TVC. 
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Figure 4. Central part of the ^^CO-observations of Dame et al. (2001). White Hnes sketch spiral arm ridges in the data. 



Figure 5. The distribution of gas in the standard gas model at 
evolutionary age 0.32Gyr. Note that the initial particle distribu- 
tion in the simulation ends at i? = 8kpc, producing the artifical 
outer cutoff of the particle disk in the plot. The position of the 
Sun at (x,y) K, (7.5 kpc, 2.7 kpc) is shown by the symbol. 



Figure 6. The terminal velocity curve (TVC) of the standard 
gas model at evolutionary time 0.32Gyr (top and bottom curves), 
and of model halo also at time 0.32Gyr (middle curves), compared 
to the HI and CO data. Model velocities have been scaled by a 
model-dependent factor g to best fit the observed terminal veloc- 
ities for 10° < \l\ < 50°. The observed ve oc ties are corrected to 

' '. The TVCs of model 



the Hipparcos LSR frame as described in § 2.4 



halo are offset by 60km/s for better readability of the diagram. 



Important features in the observed {l,v) diagram are the 
ridges of emission which indicate the location of spiral arms. 
The white lines drawn in Fig. ^ reproduce approximately 
the locations of these spiral arm ridges. These lines are 
then transformed to the Hipparcos LSR [Figure ^ was con- 
structed assuming a solar motion of \vq\ = 20kms"^ to- 
wards {l,b) = (56.2°, 22.8°)), and are then overplotted on 
most model {l,v) diagrams. These observed ridges should 
be reproduced by a good model of the MW gas flow. The 
reverse need not always be true, however, because the visi- 
bility of a spiral arm ridge in the data may depend on the 
radial distribution of gas and the geometry of the arm with 
respect to the line-of-sight. 

In the model {l,v) diagrams we also show molecular 
cloud observations of Dame et al. (1986) and Bronfman 



(1989) [symbol "x" in the p lots] , and HII region o bserv ations 
of Georgelin & Geo rgelin (1976), Downes et al. (1980), and 
Caswell & Haynes (1987) [symbol Most of our models 



do not contain a halo potential, and therefore underestimate 
the velocities for |Z| ^ 40 — 50°. Thus we omit the molecular 
cloud and HII observations in these models for \l\ > 40°. For 
the sake of clarity we have also left out clouds with less than 
from the Bronfman et al. data, as well as clouds in 
the smallest brightness bin from the Georgelin & Georgelin 
sample. For a more detailed discussion of these observations 
see Paper I. 



3 BEST-FIT MODEL FOR THE MILKY WAY 

We have investigated a number of gas flow models in the 
COBE potentials of Section |2|, for different pattern speeds, 
bar angles, and stellar spiral arm morphologies. The analy- 
sis of these model sequences is deferred to Section ^. Here 
we begin with a description of our best model for the gas 
dynamics in the MW. This best-flt model (hereafter called 
"standard model") is based on the standard four- armed 
ifi-har ~ 20° luminosity model (Bissantz & Gerhard (200' 



- Paper II), as described in Section 2, it is point-symmetric, 
the bar pattern speed is Qp = 60Gyr"\ implying corota- 
tion at _Rci ~ 3.4 kpc, and the spiral arm pattern speed is 
Qsp = 20Gyr-\ 

The gas distribution of the model is shown in Fig- 



ure ^ In this plot the Sun is at (a;, y) = (7.5, 2.7) kpc, with 
Ro ~ 8 kpc. Both inside and outside corotation there exist 
four spiral arms which are connected in a complicated way 
through the corotation region of the bar. Outside corotation 
the spiral pattern in the gas response consists of a pair of 
strong arms and a pair of weaker arms. 

To compare the model to the gas observations, we scale 
it to the observed terminal velocity curve (TVC) in the lon- 
gitude range 10° < |/| < 50°, using an LSR circular velocity 
Vo = 220 km s"'. With this scaling the CO TVC is repro- 
duced well by the model (Figure ^ , including most distinct 
"bumps" in the observations except that at / ~ 50°. Vice- 
versa, the model TVC shows an extra bump at I « —15° 
which is not seen in the data. Presumably the potential in 
this bar-disk transition region is not accurately modelled - 
there is little information in the NIR data on the mass dis- 
tribution in this region. Overall, however, the fall-off with 
longitude and even the detailed form of the inner disk Galac- 
tic TVC are represented well by the model. 

In the bulge region, for |/| < 10°, the model velocities 
are limited by the resolution of the hydrodynamic simula- 
tion. We thus do not expect to fit the observed large ter- 
minal velocities there, but several other effects may play a 
role as well. See §3.1 below and Englmaier & Gerhard (199£) 
for a more detailed discussion. Because the gas particles in 
the numerical simulations flow inwards, however, the veloc- 
ity structure in the inner 1 — 1.5 kpc is not of significant 
relevance for the gas flow and, particularly, for the shock 
structure in the main spiral arms well outside this region. 
Only the velocities of the pair of inner arms passing the mi- 
nor axis of the bar laterally at « 1 kpc could be somewhat 
affected. Since the main aim here is to investigate the large- 
scale gas flow and spiral arm morphology, we have therefore 
not attempted a detailed fit to the inner bulge terminal ve- 
locities. 

The standard model does not contain a dark matter 
halo. For the assumed LSR velocity of Vo = 220 km s^^, 
its TVC falls below the observed TVC at |/| > 40° - 50°. 
At \l\ — 90° the deficit in the model's terminal velocity is 
~ 25km/s. This can be corrected by a adding a halo poten- 
tial that generates a rotation velocity vdm ~ 120km/s near 
the Sun. We thus constructed a new gas model in a poten- 
tial which includes a suitable quasi-isothermal dark matter 
halo, which has circular velocity at infinity Vi„t = 220km/s 
and core radius a = 10.7 kpc (there is substantial freedom in 
these parameter values) . The TVC of this model halo is also 
shown in Fig. ^ its best-fitting scaling factor ^ and mass-to- 
light ratio is only slightly different from that of the standard 
model. Model halo fits the shape of the observed TVC out 
to _Ro, but is not quite as good a match to the bumps in 
the TVC as the standard model. This suggests that the as- 
sumed halo model is oversimplified. Note the fact that the 
constant M/L standard model, a maximum disk model by 
construction, reproduces the observed terminal velocities in- 
side galactocentric radius R « 5 kpc and still accounts for 
most of the circular velocity near the Sun. 

The {l,v) diagram for the standard model is shown in 
Figure ^. The dense ridges in this diagram show the loca- 
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Figure 7. (/, v) diagram for the standard gas model at evolutionary age 0.32Gyr, for LSR-velocity Vo 
are scaled with the factor ^ determined from Fig. ^ 



: 220km/s. The model velocities 



Figure 8. Grey-scale Zu-plot for the standard gas model with a dark halo component included in the potential (model halo), at 
evolutionary age 0.32Gyr. Particles in regions of low gas surface density {< 1% of the maximum surface density) are suppressed to 
enhance the contrast, mimicking also the bias of the observed distribution of molecular gas and HII regions towards higher densities. 
Model velocities are scaled such as to fit best the observed terminal velocities in 10° < \l\ < 50° (see Fig. H); the LSR-velocity is 220km/s. 
For comparis on wi th observations, t he sp iral arm ridge lines from Fig. ^are overplotted, as are the data for gian t molecular cloud s from 
Dame et al. ( |l986[ ) and Bronf man (L989) ("x" symbols), and for HII regions from Georgelin & Georgelin ( |l976| ), Downes et al. ( |l98C| ), 
and Caswell & Havncs ([igSTj ) symbols); see ^.4| . The short vertical lines mark the observed spiral arm tangent directions from 

Englmaier & Gerhard (1999). 



Figure 9. Contours of constant radial velocity for gas clouds 
in model halo, as seen by an observer moving with the velocity 
of the LSR. Contours are spaced by lOkm/s. Dashed contours 
indicate negative radial velocities, full contours positive radial 
velocities. Ticks along the full line give distances from the LSR 
(right end), in kpc along the line-of-sight through the Galactic 
Center. Longitudes as seen from the LSR are indicated on the 
margin of the plot. 



Figure 10. A sample of closed xi and X2 orbits in the bar frame, 
in the ipbar = 20° potential of our standard gas model with pat- 
tern speeds Qp = 60Gyr~^, Q^p = 20Gyr~^. Note the convergence 
of the outer X2-orbits on the major axis of the bar. 



tions of the spiral arms. Figure 11 of Paper II shows the 
correspondence between the locations of spiral arms in the 
Galactic plane and in the (/, v) diagram. The model spiral 
arm ridges in Fig. ^ generally coincide very well with the 
observed spiral arm ridges. Less good is the correspondence 
for the 3kpc-arm (again in the bar-disk transition region), 
which is at too small negative velocities in the model com- 
pared to the observations for I > —5°. 

For model halo the Iv-plot is shown in Figure ^ at en- 
hanced contrast to emphasize the spiral arms (see 



2.4) 



Compared to the standard model {I, v) diagram (without 
halo), significant differences are in the outer spiral tangents 
near I « ±50° which are relocated by a few degrees, and 
in the terminal velocities at |/| > ±50°. This is a general 
result: For a number of models discussed later in this pa- 
per we have added a suitable halo potential and computed a 
new gas model, sometimes additionally changing the extent 
of the initial gas disk in the simulation from the standard 
8 kpc to 10 kpc. In all these cases the only significant change 
has been that the outer tangent points moved outwards in 
longitude by |AZ| < 5°. Thus we do not include a halo po- 
tential in the remaining models discussed below. 

Also shown in Fig. ^ are the observed spiral arm ridge 



lines, tangent points, and tracers from {2.4. The comparison 



with the grey scale plot for the model shows that model halo 
gives a very good description of the gas kinematics in the 
disk outside the bar. For many features deviations are less 
than ~ lOkms"^. 

We end this section with showing in Figure ^ a map of 
the radial velocities of gas clouds with respect to the LSR for 
model halo. This map allows one to assess the likely errors 
made in determining kinematic distances from cloud radial 
velocities by assuming a circular orbit model. 



Figure 11. The same orbits as shown in Fig. now displayed 
in an {I, v) diagram, using i/Jbar = 20° also for the projection. The 
orbit with the highest peak velocity, at ( ~ 2°, is the cusped orbit. 
Inside this cusped orbit two xi-orbits with self-intersecting loops 
are plotted (cf. Fig [lo[). Slight oscillations in the (l,v) traces of 
the outer xi orbits betray the time-dependence of the potential. 
All orbit velocities are scaled by the same factor t; as in the TVC 
of the standard model in Fig. H. Also shown in the figure are the 
observed terminal velocities in the bulge region. The cusped orbit 
and the other xi -orbits at larger galactocentric radius represent 
the observed terminal velocities well. The velocities of X2-orbits 
peak at \v\ ~ 85km/s in the plot. 



3.1 Orbits and gas flow in the bulge region 

In this section we consider in more detail the gas flow and 
TVC in the central \l\ < 10°. The discussion of the large- 
scale morphology and pattern speed is continued in 

We begin with an analysis of closed orbits in the stan- 
dard potential, including a c entral cusp in the mass model 
as described in Section 2.2, and assuming the same val- 
ues for the pattern speeds as in the standard gas model. 
The closed orbits are found with a simple shooting algo- 
rithm. Despite the intrinsic time-dependence of the poten- 
tial in the bar frame, the closed orbits in the inner kpc re- 
main essentially unperturbed. Figure ^ shows closed X2- 
and xi-orbits around the so-called "cusped" orbit, the xi- 
orbit whose turning points on the bar's major axis have a 
cusp shape. Closer to the galactic centre the more tightly 
bound xi-orbits become self-intersecting, and the X2-orbit 
family of stable orbits elongated perpendicular to the bar 
appears. Note that in this potential the outermost X2 orbits 
are (nearly) converging on the major axis of the bar, imply- 
ing that gas clouds on these orbits would collide. This limits 
the radial extent of accessible X2-orbits for gas clouds in a 
hydrodynamic flow. 

Between the cusped xi orbit and the first non- 
intersecting X2 orbit there is a region in which no closed 
orbits suitable for gas fiow exist (Fig. ^o|). Because of this, 
infiowing gas has to quickly pas s this region: the mecha- 
nism described by Binney et al. (1991), where gas moving 



in from the last ii-orbits collides with gas on the outer- 
most a;2-orbits, producing a spray that then forms an off-axis 
shock by hitting the far side of the a;i-orbits, cannot work 
as well in the potential here because of the lack of suit- 
able outer a;2-orbits. Instead, the main place of dissipation 
of kinetic energy is likely to be the self-crossing loops of the 
ii-orbits inside the cusped orbit, from where the gas moves 
inwards to hit the a;2-disk further in. This may explain why 
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there is a mostly gas-free gap in the hydrodynamic simula- 
tions between the last non-intersecting xi orbit and the first 
accessible X2 orbit, without any clear off-axis shocks like 
those found in other barred galaxy models (e.g. Athanas- 
soula ( |f992[ )). 

How would gas clouds following the closed orbits in our 
standard model potential compare to the terminal velocity 
observations? Figure ^ plotted for bar angle 20°, shows 
that the terminal velocities of the closed Xi-orbits repro- 
duce the observations surprisingly well. The cusped orbit in 
Fig. |l^ is not only at the same longitude as the maximum 
in the observed TVC (at positive / where we have data), but 
also appears to account for the decline of the observed termi- 
nal velocities at lower I. This is consistent with a gas- free gap 
between the cusped orbit and the first acceptable X2-orbit, 
due to which no strong leading shocks form like those ob- 
served in other barred galaxies. Episodic infall of gas clouds 
within the gap region may nonetheless form transi ent sh ocks 
similar to those observed by Hiittemeister et al. (1998). 



from NIR surface bri ghtne ss photometry of early type spi- 
rals, Moriondo et al. (1998) obtain a mean value and disper- 



The closed orbits thus reproduce the observed terminal 
velocities, but why is there no gas at these velocities in the 
simulations? Possible explanations are as follows: First, the 
resolution length of our SPH code, which is dominated by 
the smoothing length of the SPH particles, may be too large 
to follow the strongly elongated xi-orbits near their cusped 
ends. Also, in the low density region further in, the relative 
velocities of neighbouring particles are large, so that the cor- 
responding viscosity may lead to fast infall of SPH particles 
to the centre. This would depopulate the inner xi orbits, 
where the largest terminal velocities are expected. Indeed, 
the terminal velocity curve attains larger peak velocities in 
high-resolution models with some lO'^ particles (the stan- 
dard gas model has ~ 6 ■ 10^ particles), for example the 
maximum is at « 222km/s on the I > 0° side, compared to 
« 210km/s in the standard gas model. 

Second, the detailed orbit shapes in the inner kpc of 
the MW depend on the fine-structure of the potential there, 
which is not accurately known, both because of the limited 
resolution of the underlying NIR data and luminosity model, 
and the difficulty of the deprojection in this region. The true 
xi orbits could easily be somewhat less cuspy near their 
ends, or the a;2-orbits less converging, making them more 
eeisily populated with gas clouds. 

Third, the hydrodynamic flow in the central MW could 
be genuinely slower than suggested by the closed orbits 
which it approximately follows. Increasing the model ter- 
minal velocities in the bulge region to match the observa- 
tions would then require a somewhat larger bulge mass-to- 
light ratio than the value determined from the TVC fit at 
10° < |Z| < 50°. This might be plausible if the disk popula- 
tion is somewhat younger than that of the bulge. It is un- 
likely that the M/L ratio is significantly larger in the entire 
central kpc of the MW because the peak velocitiy of gas on 
X2-orbits in our standard gas model is 85km/s, which com- 
pares well to the observed < 80km/s for CS-cloud cores at 
|/| < 0.7°, where the projected model X2-orbits are located. 
However, the potential near 7? ~ 1 kpc, which is most rele- 
vant for the xi orbits, is more sensitive to the upper bulge 
component than that in the inner 100 pc where the X2 orbits 
are. The required variations of the bulge M/L to give termi- 
nal velocities of > 250km/s would seem consistent with the 
spread of K band mass-to-light ratios for other bulges. E.g., 



sion (0.6 ± 0.2) (M© / Lq)k [the velocity scale of the standard 
gas model corresponds to {M/L)k ~ 0.6{Mq/Lq)k]- 

In summary, there are uncertainties in modelling the 
gas velocities in the inner few hundred parsec of the bulge. 
We will not pursue this further here, because as already 
discussed in the last subsection, this is not important for 
the gas flow in the main spiral arms further out. Rather, we 
now turn to the determination of the pattern speeds in the 
Milky Way. 



4 THE GAS FLOW IN MODELS WITH 
SEPARATE BAR AND SPIRAL ARM 
PATTERN SPEEDS 

Figure ^ shows a time sequence for the gas distribution 
in our point-symmetric standard model with pattern speeds 
flp = 60Gyr"^, fJgp — 20Gyr"\ The inner arms emanat- 
ing from the ends of the bar corotate with the bar pattern 
speed; they are clearly driven by the bar. These inner arms 
are connected to the outer spiral arms by a time-dependent 
transition region near bar corotation. Here a lateral arm 
from the distant end of the bar as seen from the Sun merges 
with the other inner arm from the nearer end in some frames 
(4-5), but not in others (1-2). In the latter case, the lateral 
arm continues into the outer spiral arm passing close to the 
Sun, in the former both arms join into a weaker outer arm 
staying well inside the Sun. The outer arms themselves move 
with respect to the bar frame. However, they do not rotate 
in the plot steadily around the model's centre, as one might 
have expected if they were driven by the different pattern 
speed of the spiral arm potential, but apparently exhibit 
complicated back-and-forth oscillations, with respect to each 
other and with respect to the bar frame, and some arms 
merge and bifurcate at certain times; compare the vicinity 
of {x,y) — ( — 5, —3) in the different panels. 

The outer spiral structure evolves in such a complicated 
way because both the spiral arm potential and the bar simul- 
taneously force the gas distribution with different pattern 
speeds. To investigate this further, we have computed a gas 
model similar to the standard model, but with spiral struc- 
ture removed from the gravitational potential, so that the 
non-axisymmetric component of the potential is solely due 
to the bar. The resulting model has four spiral arms inside 
the bar corotation radius, and two symmetric outer spiral 
arms, all stationary in the bar frame. In the model with a 
driving spiral structure potential, the four outer arms can 
be regarded as a superposition of one component generated 
by the bar perturbation, which is at nearly constant posi- 
tion in the bar frame, and a second component driven by 
the spiral structure potential perturbation. The former is 
the stronger pair of arms in Fig. |l^, one of which passes just 
inside the Sun symbol in the flgure. The second component, 
visible as the weaker pair of arms between the bar-driven 
arms in Fig. ^| can be seen to fall behind the bar-driven 
arms through the sequence of frames in Fig. |l| {flpyflsp). 
At certain times the second component is seen to branch off 
the bar-driven arms (frame 1). Thereafter, the spiral-driven 
arm appears to fall behind, move inwards, until it finally col- 
lides with the opposite bar-driven arm (frames 2 through 6 
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Model 


Potential 


f2p 






Remarks 




standard 


UUVjry 1 


1 


1 Q F,^^^.rr^^ 

ly.u'^yi 


otayiidajro. moQel 


40 


standard 


61.4Gyr 


-1 


40.8Gyr-i 




60 


standard 


61.4Gyr 


-1 


61.4Gyr-i 




halo 


standard 




1 


in fir^Trf— 1 
iy.DLryr 


btandard model 












ivi/^lii/^ivirr nal/^ 


barSO 


standard 


50Gyr- 


1 


20Gyr-i 




bar70 


standard 


70Gyr- 


1 


20Gyr-i 




opGn2 




60Gyr" 


1 


20Gyr^^ 






to standard four-armed model 










2spi 


two-armed, pitch angle half 


60Gyr- 


1 


20Gyr-i 






the value in the standard four-armed model 










mix 


similar to four-armed model, but Sag-Car arm 


eOGyr- 


1 


20Gyr-i 






and counter-arm are reduced in amplitude 










lHJfXL Ills 


Q'f"Civi/^QT*/H ■^st'tH^ ct^ttqI t\£^v4" ii"¥'r\Q4"ir^Ti cnrT"t"r*n^/H /^tt 

t) LdliLlCLl U. VViljll a IJll di JJCl L Lll UdLiUll O W i LUlltJU. KJll 


60Gyr" 


1 


20Gyr^^ 




incllO 


four-armed model, bar angle 10° 


60Gyr- 


1 


20Gyr-i 




incll5 


four-armed model, bar angle 15° 


60Gyr- 


1 


20Gyr-i 




incl25 


four-armed model, bar angle 25° 


60Gyr- 


1 


20Gyr-i 




incl30 


four-armed model, bar angle 30° 


60Gyr- 


1 


20Gyr-i 




strongarms 


standard 


60Gyr- 


1 


20Gyr-i 


m > 2-multipoles of density 












outside 3.5 kpc multiplied by 1.5 


tumblingbar 




60Gyr- 


1 


20Gyr-i 


standard potential, but the centre 



of the bar perturbation is not in MW centre 



Table 1. The gas models discussed in this paper. 



Figure 12. The gas distribution of the standard model 20 for a sequence of evolutionary times and corresponding phase differences 
between the bar and spiral components of the potential: (0.296Gyr, 318°), (0.304Gyr, 337°), (0.312Gyr, 356°) (upper row from left to 
right), and (0.320Gyr, 14°), (0.328Gyr, 33°), and (0.336Gyr, 51°) (lower row from left to right). The long axis of the bar is aligned with 
the X-axis in all panels. Note the evolution in the connecting region between the inner and outer arms. 



~ frame 1 through frame 3) . In the course of this evolution, 
the detailed morphology of the transition region around bar 
corotation changes. At most times, the spiral-driven arms 
are connected to the inner lateral arms (frames 2-5), at oth- 
ers the latter connect to the bar-driven arms (frame 1). At 
certain times, an arm may look fragmented, and its tangent 
point may split in longitude. 

Thus when neither the bar nor the spiral structure com- 
ponent dominates the non-axisymmetric potential and gas 
flow in the disk, as in Fig. it is difficult to deduce from a 
single snapshot of the arm morphology that the system sup- 
ports a second, independent pattern speed. This may explain 
why it is difficult to say from the observed arm morphologies 
in barred spiral galaxies whether the spiral ar ms ar e driven 
by the bar or not (see Sellwood & Wilkinson (1993)). 



4.1 Pattern speeds in the Milky Way 

We now investigate Galactic models with different combina- 
tions of the bar pattern speed Qp and the spiral arm pattern 
speed flap- Particularly for the spiral arms the assumption 
of a constant pattern speed is probably still idealized, but 
represents a first step towards understanding realistic cases. 
In these models, the separation of bar and spiral arm compo- 
nents in the potential is based on the density distribution. 
We assume that the m > 2-multipoles of the density at 



galactocentric radii r < Tcut = 3.5 kpc (bar) rotate with Q,p, 
and outside of rcut (spirals) with ilap- The specific value of 
Tcut was chosen because it corresponds approximately to the 
end of the bar in the reference luminosity model of Paper 
II, and is also equal to the inner radius of the spiral pattern 
there. All models described in this Section are based on this 



luminosity model (cf. §2.1) 



Multiple pattern speeds {flp flap) imply a genuinely 
time- dependent potential in the frame corotating with the 
bar. For each model we have therefore investigated a se- 
quence of snapshots at different evolutionary ages, separated 
by A0 ~ 23° in phase difference between the two compo- 
nents in the potential. This corresponds to steps in the evo- 
lutionary age of the model of « 0.01 — 0.02Gyr, depend- 
ing on the combination of pattern speeds. A finer analy- 
sis of our standard gas model in steps of A(p « 9° showed 
that no significant features in the models are missed with 
A(j} « 23°. Because the potentials used in this section are 
point-symmetric, we need to cover only a range of 180° in 
bar-spiral phase difference. All analysed snapshots have evo- 
lutionary ages at or around 0.30Gyr, the time after which 
the gas flow in the similar single pattern speed models of 
Paper I had become approximately quasi-stationary. 

We begin with the three models 60, 40 and 20 with 



different spiral arm pattern speeds Qs 



61.4,40 and 



19.6Gyr . The bar pattern speed is set to 61.4Gyr 
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Figure 13. Gas distribution in models 60 (left), 20 (middle), and 40 (right). Overplotted black dots indicate particles highlighted in the 
(/,!)) diagrams in Fig. ha. The position of the Sun at {x,y) (7.5 kpc, 2.7kpc) is shown by the symbol. 



Figure 14. Terminal velocity curves for the best evolutionary 
ages of the single and double pattern speed models 60, 20, and 
40. The model TVCs are scaled by different factors ^, determined 
respectively from the best fit to the observations for 10° < |i| < 
50°. For clarity, the TVCs are separated by 40km/s from each 
other in velocity. 



(60Gyr~^ for model 20), based on the results of Paper I. 
Model 60 is the single pattern speed model; the two other 
values for Qap are motivated as follows. Firstly, Amaral & 

20Gyr-\ 



Lepine (1997) preferred the lower value D,sp 
based on the positions and ages of open clusters, and from 
compa ring their model with other data. Fernandez et al. 
(2001) obtained a somewhat higher ilsp ~ 30km/s/kpc from 
Hipparcos data for OB stars and Cepheids, supplemented 
with radial velocities or distances, respectively. Thus we se- 
lected flap ~ 40Gyr"^ as an additional intermediate value 
between i}sp = 20Gyr^^ and the single pattern speed model. 
Secondly, for Qsp = 60, 20, 40 Gyr~^, the bar corotation 
radius coincides approximately with the corotation, inner 
Lindblad, and inner ultra- harmonic (1:4) resonance of the 
spiral pattern, respectively (see Fig. 

What are the characteristic differences between these 
models? In Figure |l3| we compare their spatial gas distribu- 
tions. All three models have four inner spiral arms inside bar 
corotation (at ~ 3.4 kpc), and four outer spiral arms outside 
of a transition region beyond the end of the bar. Because we 
have always selected that snapshot of a model which best 
reproduces the observed {l,v) diagram (Fig. ^), the spiral 
arm tangent points are fixed relative to the position of the 
Sun. Thus in the figure the outer arms of all models appear 
approximately at the same positions. 

However, near the bar corotation radius the models dif- 
fer significantly. In the single pattern speed model 60, where 
the spiral structure corotation radius coincides with that of 
the bar, the gas in the corotation region moves with near- 
sonic velocities relative to the pattern. Consequently shocks 
are weak or nonexistent in this region, and the spiral arms 
in the gas response dissolve there. On the other hand, for 
model 40 with flap — 40Gyr~^, the corotation radius of the 
spiral structure is at ~ 5 kpc. Again, the gas distribution 
shows gaps in the spiral arms at this location - one pair 
of arms nearly vanishes there, the other pair weakens - but 
there is no such gap in the spiral arms near bar corota- 
tion. Rather, the connection between the inner and outer 
spiral arms is dynamic, due to the different pattern speeds. 
Finally, in model 20 the spiral structure corotation radius 
occurs beyond the solar orbit. In this model we see no clear 
gaps in the arms. However, near the bar corotation radius 

~ 3.4 kpc, which coincides with the inner Lindblad res- 
onance of the spiral pattern, the arms weaken, and the tran- 
sition region appears to be more complicated than in model 
40. 

The comparison with the Galactic TVC is shown in 
Figure ^ All three models reproduce the observed data 
quite well, but the slope in the region 10° < \l\ < 50° is 
fit significantly better by the multiple pattern speed models 



than by the single pattern speed model. For |/| > 40° — 50° 
the model TVCs start to fall below the observed TVC; as 
we have discussed in for the assumed Vo = 220kms~^ 
the dark halo starts to contribute to the observed velocities 
there, according to our NIR-based models. 

Figure [Tsl shows {I, v) diagrams for the three models. To 
facilitate their interpretation, we have overplotted a number 
of observed features, and have highlighted particles belong- 
ing to specific arm features both in this figure and in the 
density plots of Fig. ^ labelled by "A" and "B". Feature 
"A" corresponds to the 3kpc-arm in the observed CO {l,v) 
diagram (Fig. ^ , feature "B" to its symmetric counter-arm 
in the models (see Fig. |l3[ in model 20 the neighbouring 



arm "A' " is included). 



Particularly important is the existence of certain re- 
gions in the {l,v) diagram, e.g., next to the 3kpc arm (fea- 
ture "A"), where hardly any gas is found in the multiple 
pattern speed models. Such voids, designated by "V" in the 
plots, arise because the gas is aligned morphologically and 
kinematically with spiral arms there, which appear as well- 
defined ridges with adjacent voids in the {I, v) diagram. From 
Fig. ^ we see that similar voids are also visible in the ob- 
served ^^CO {l,v) diagram, specifically next to features "A" 
and "B" in the bar corotation region. This suggests strongly 
that in the MW the gaseous arms go through bar corotation. 
On the contrary, in the single pattern speed model 60, the 
arms dissolve in the bar corotation region, the 3kpc-arm is 
thus imcomplete, and the gas is spread out approximately 
evenly over the corresponding parts of the {I, v) diagram. 
Models with a separate second pattern speed for the spiral 
arms in the MW are therefore preferred over single pattern 
speed models. Models with a growing bar amplitude also 
suppor t spir al arms in the corotation region (Thielheim & 
WolflF (|l98^)); however, it is likely that when self- gravity is 
included and the amplitude becomes non-linear, the grow- 
ing spiral pattern will again develop an independent pattern 
speed. 

Overall, the models in Fig. |l^ provide a good match 
to the observed CO (Z, v) diagram. The main features that 
are reasonably well represented are: the arm tangent at 
I — 30°; the observed I — 25° tangent, although at smaller 
I ~ 20° in models 20 and 4O', the morphology of this 
arm; the location of the main spiral arm leading to the 
/ = 50° tangent; the morphology of the ridges and voids at 
— 10° > I > —25° (apart from model 60); the arm morphol- 
ogy around {I, v) — (20°, 60kms^^). The main weaknesses of 
the models, if we disregard the poorly resolved bulge region, 
are: the arm in the models which, returning from the I = 20° 
tangent, crosses the region 5° < Z < 10° at u > 50km/s and 
does not have a counterpart in the data, except perhaps if 
shifted to lower velocities; the missing tangent at Z = —30° 
in models 60 and 20; the displacement of the 3kpc-arm to- 
wards lower velocities. 

In the model 20, the spiral arm corresponding to the 
Centaurus tangent at Z ~ —51° is strongest when one of the 
spiral-driven arms coincides with the bar driven-arm (e.g., 
frames 3-4 in Fig. [l^ ). Shortly before and after this evolu- 
tionary time this arm looks fragmented in the model, and 
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Figure 15. {l,v) diagrams for tlie best evolutionary times of models 60 (at age O.SOGyr, top), 20 (at 0.32Gyr, middle), and 40 (also at 
0.32Gyr, bottom). The model velocities are scaled by the same factor § as the respective TVC in Fig. the LSR-velocity is 220km/s. 
Particles in low surface density regions are suppressed to enhance the contrast as described in the caption of Fig. ^. For comparison with 
observations, spiral arm ridge lines from Fig. M, positions of giant mo lecula r clouds ("x" symbols) and HII regions ("+" symbols) from 
Fig. jsj and spiral arm tangent directions from Englmaier &i Gerhard (1999) are overplotted. 



Figure 16. {I, v) diagrams for the best evolutionary time 0.307Gyr of model barSO (top, Qp 
fear70 (bottom, Qp = 70Gyr~^), with model velocities scaled as described in the caption of Fig 
Figs. |and|l5[ 



50Gyr-i), and for 0.321Gyr of model 
pi Overplotted data are the same as in 



its tangent point appears split in longitude. Indeed, there 
are observational indications that the Centaurus tangent is 
split into two parts at I « -50° and « -55° - (-58°) (cf. 
Table 1 of Paper I). 

The displacement of the 3kpc-arm in the model and the 
apparent absence of its counter-arm in the data might indi- 
cate a non-point-symmetric mass distribution i n th is region. 
This possibility is discussed further in Section 5.4, Also, in 



models where the spiral pattern speed differs from the bar 
pattern speed, we expect also the mass distribution in the 
transition region to be generally time-dependent. Thus in 
this region the mass distribution based on the NIR data can 
represent only one snapshot in time. In addition, there is not 
much information in the NIR data on the spiral arm heads 
in this crucial transition region, to constrain the luminosity 
model of Paper II. Hence the potential in this region is likely 
to be at best approximately correct. 

In summary, the multiple pattern speed models repro- 
duce the observed features in the bar corotation region bet- 
ter than the single pattern speed model, in particular the 
regions void of gas in the {I, v) diagramnear the 3 kpc arm. 
Comparing the Iv-plots for models 20 and 40, we have found 
a slight preference for model 20 from the positions of the spi- 
ral arms; however, the differences between these two mod- 
els are too small to determine the spiral arm pattern speed 
rehably. We take the model with flsp — 20Gyr"^ as our 
standard model, and use this pattern speed in the following. 

We now proceed to determine the best bar pattern 
speed Qp. To this end we consider two models barSO and 
harlO which have fip = 50Gyr~^ and 70Gyr~^, respectively. 
For both models, we set Q.sp = 20Gyr"^ ; since the differences 
between models 20 and 40 are small, the precise choice of 
S7sp should not matter. Figure ^ shows {l,v) diagrams for 
the respective best evolutionary times of models barSO and 
barlO. For both models the fit of the spiral ridges in the 
inner disk region of the (/,«) diagram is significantly worse 
than for the models with Sip ~ 60Gyr^^. In model barSO the 
positions of the spiral arm ridges and tangents are worse 
(e.g., one of the I > tangents is missing), and the I < 
TVC and the 3kpc-arm are particularly badly represented. 
In model barlO the spiral arm ridge with tangent position at 
I « 25° is at significantly too small longitudes, and the same 
is true for the Centaurus tangent near / ~ —51°, which cor- 
responds to one of the bar-driven outer arms (see discussion 
above). It is noteworthy that the position and shape of the 
Centaurus tangent are sensitive to S7p despite a galactocen- 
tric radius of this tangent point of more than 6 kpc. Note 
also that a pattern speed of 70Gyr"^ would put corotation 
inside the end of the NIR bar (see §^ and Paper II). 

We conclude that the best value for the bar pattern 



Figure 17. TVCs at the respective best evolutionary time, for a 
sequence of models with bar angles ySbar = 10° (model incllO, 
at age 0.32Gyr), i/^bar = 15° (model incll5, at age 0.32Gyr), 
ipbar = 25° (model incl25, at age 0.31Gyr), and ipbap = 30° (model 
inclSO, at age 0.30Gyr). For comparison, the standard model 20 
for bar angle i/^tar = 20° is included in the figure. In all models the 
pattern speeds are Qp = 60Gyr^^ and Qsp = 20Gyr^^, and model 
velocities have been scaled by factors 5, determined for each model 
by fitting the observed terminal velocities for 10° < \l\ < 50°. For 
clarity, these TVCs have been offset in steps of 40km/s, with 
model incll5 plottet at the correct velocities. 



Figure 19. TVCs at the respective best evolutionary time, for 
gas models forced by different types of spiral structure in the 
mass density. Model open2 is shown at age O.SOGyr, model 2spi 
at age 0.30Gyr, and model mix at age 0.31Gyr, model strongarms 
at age 0.32Gyr. Also shown is the standard model 20, and model 
noarms without any massive spiral arms. For all models, the pat- 
tern speeds are Qp = 60Gyr~^ and Qsp = 20Gyr~^, and the ve- 
locity scale is fixed by fitting to the observed terminal velocities 
for 10° < |i| < 50°. For clarity, the model TVCs are offset in steps 
of 40km/s, with model open2 plotted at the correct velocities. 



speed in the Milky Way is Sip = (60 ± 5)Gyr"\ For the 
preferred scaling of the employed NIR bar model this corre- 
sponds to bar corotation at 3.4 ±0.3 kpc, equal to the length 
of the bar within the uncertainties. I.e., the Milky Way bar 
is a fast bar. 



5 BAR ORIENTATION AND INFLUENCE OF 
STELLAR SPIRAL ARMS 

5.1 Bar angle 

In Paper II, luminosity models were generated from 
the COBE/DIRBE L-band data for bar angles (ySbar = 
10°, 15°, 20°, 25°, 30° and 44° in the same way as for the 
standard tpbai ~ 20° model. From considering the photo- 
metric residuals and the line-of-sight distributions of clump 
giant stars in the bulge, a preferred range for the bar angle 
15° < <^bar < 30° was found. Here we obtain independent 
constraints on (p-tar from corresponding gas fiow models. 

For each value of the bar angle, the luminosity model 
was converted to a mass model assuming constant M/L. As 
described above, the non-axisymmetric part of the potential 
was split into bar and spiral arm components, and gas mod- 
els were computed with pattern speeds = 60Gyr~^ and 
Slap = 20Gyr~\ respectively. In Figures h/a and |l8| we show 
the TVCs and {I, v) diagrams of models incllO, incllS, incl25 
and incl30 with bar angles ip-h^i = 10° , 15° , 25° , 30° at their 
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Figure 18. (l,v) diagrams for tlie gas flow models with bar angles 10° (top left), 15° (top right), 25° (bottom left) and 30° (bottom 
right), for the same evolutionary times as in Fig. Il7|. 



respective best evolutionary times, as judged by comparing 
their {l,v) diagrams with Fig. ^. 

Model incl25 reproduces the terminal velocity observa- 
tions with similar quality as the standard gas model 20, but 
model inclSO is clearly inferior - the terminal velocities are 
too large at positive longitudes and too small in modulus at 
negative longitudes. Obviously, this cannot be corrected for 
by a change of the velocity scaling factor ^. The TVCs of 
models mcllO and incllS axe also inferior fits to the observed 
terminal velocities, albeit not as bad as model incl30. 

The v) diagrams of these models depend strongly on 
the assumed bar angle. Models mcllO, incllS and incl30 do 
not reproduce the observed spiral arm ridges. Models incllO 
and incllS also do not have a well-define 3 kpc arm at posi- 
tive longitudes, and model incl30 shows very little arm struc- 
ture at all in the region around bar corotation. Model incl25 
reproduces the observations for a large range of longitudes 
similarly well as the standard model, but the positions of the 
spiral arm ridges at —30° ^ I ^ —10° are not matched well 
and the non-circular velocities in the 3 kpc arm are smaller 
than in the standard (pbar = 20° model. From the {l,v) di- 
agrams model 20 is best, but we consider model 25 as still 
satisfactory. 

One may ask, how much of the difference between these 
models is due to the different shapes of their underlying lu- 
minosity distributions and gravitational potentials, and how 
much of it is due simply to the different viewing geome- 
tries with respect to the gas ffow? To answer this question, 
we have constructed the {l,v) diagram of the standard gas 
model 20 seen from a viewing angle of (^bai ~ 30°. The mor- 
phology of the arms in this {I, v) diagram differs little from 
that in the original y>bar ~ 20° {l,v) diagram because of 
the tightly wound pattern. Thus the differences between the 
{I, v) diagrams in Fig. ^ must be caused mainly by genuine 
changes in the gas flow, originating from the different mass 
distributions corresponding to the COBE data for different 



5.2 Spiral arm models 

The MW very probably has four spir al arrn s in the gas dis- 
tribution (e.g., Englmaier & Gerhard ( 1999 )). It is not clear, 
however, whether these also correspond to four stellar spiral 
arms, because some of the tangent points are not clearly seen 
in the near-IR light (see th e disc ussion in the Introduction 
and in Drimmel & Spergel (2001)). Here we investigate the 
{l,v) diagrams that result when different stellar spiral arm 
patterns drive the gas flow, and flnd evidence for a four- 
armed spiral pattern also in the distribution of luminous 
mass. Speciflcally, we have studied three models: 



Open2: A two-armed model with the same spiral arm pitch 
angle 13.8° as in the standard (four-armed) mass model, but 
where the Sag-Car arm and its counter-arm are removed 
from the model. Such a model can be justifled by the fact 
that the Sag-Car arm is hardly visible in the NIR. 

2spi: Another two-armed model, but with approximately 
half the pitch angle of the standard four-armed model. This 



model reproduces approximately the same tangent point po- 
sitions on the sky as the standard model. 

Mix: Similar to the standard mass model, but the Sag-Car 
arm and its counter-arm are given only 40% of the peak 
density amplitude of the other two arms. This is based on 
the result of Drimmel & Spergel ( 200l| ), who found that 
they had to reduce the amplitude of the Sag-Car arm in 
their best-fit four-armed model for the COBE/DIRBE J and 
K-band NIR maps. 

Luminous mass models with these patterns were de- 
rived using the algorithm described in Paper II, all for bar 
angle (p-t^i = 20° , by incorporating the respective spiral arm 
model in both the parametric initial model and in a penalty 
term for the non-parametric deprojection. As in the stan- 
dard mass model the non-parametric algorithm changes the 
spiral arms somewhat, but the overall pattern stays intact. 
Because the NIR data constrain only the arm tangent points 
(Paper II), two-armed and four-armed models which repro- 
duce the tangent point data fit the photometry with similar 
quality. 

The TVCs of gas fiow models computed in the corre- 
sponding gravitational potentials are shown in Figure |l9[ 
The overall slope with I of all these model TVCs is similar, 
because this is dominated by the monopole term in the mass 
distribution. However, all three models open2, 2spi and mix 
do not fit the wavy structure of the TVC data as well as 
the standard model 20 (Fig. with model mix the best 
among the three. For comparison we also include in Fig. |l^ 
the TVC of a model noarms, whose gas fiow was determined 
in a potential that includes only the perturbation from the 
bar, but not that from the spiral arms. This model has only 
two outer spiral arms in the gas distribution. 

{I, v) diagrams of these models are shown in Figure |20| 
Models open2 and 2spi compare poorly to the observed CO 
(l, v) diagram: the up-turning arm at / ~ 25° is missing, and 
the spiral structure is generally wrong for —15° < I < —35°. 
In both these models the fit to the data is bad because the 
gas distribution of the model is only two-armed. We have 
checked that this is not a result of choosing the wrong bar 
and spiral arm pattern speeds, by computing gas models in 
the two-armed potential of model open2 for the additional 
combinations of (fip = 50Gyr~^,n5p = 20Gyr~^) and Q,p — 
Qsp ~ 60Gyr"^ (single pattern speed). In both cases, the 
global, two-armed morphology of the gas flow is the same as 
in model open2. Model mix is the best of the three models 
in this section. However, the envelope of its Iv-plot shows 
stronger bumps than the standard model, in particular near 
—30° ... — 15°. The gravitional potential of model mix is 
quite similar to our standard potential, and therefore the 
similarity of the gas flows is expected. From these tests we 
conclude that a four-armed spira l arm potential is pre ferred 
(see also Englmaier & Gerhard ([l999|) and Fux (|l999|)). 



5.3 Stronger spiral arms 

The standard gas model fits well most spiral arms in the 
observations. However, the famous 3kpc-arm, a prominent 
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Figure 20. {I, v) diagrams of the models with TVCs shown in Fig. |l^, with their respective velocity factors ^. From top to bottom: 
model open2, model 2spi, and model mix. 

Figure 21. Lv-plot for model strongarms that has strong spiral arms in the potential, at its best evolutionary age 0.32Gyr. Model 
velocities have been scaled by an appropriate factor ^, as described in the caption of Fig. 



feature which extends from {I « 10°,?; = Okm/s), through 
(l = 0°,v^ -50km/s), to (/ ^ -22°, w « -120km/s) is still 
displaced towards lower non-circular velocities. Because the 
underlying mass model already contains spiral arms, assign- 
ing gravitating mass to the gas particles does not lead to an 
improvement (as a corresponding model confirmed). 

However, compared to near-IR-observations of other 
spiral galaxies (cf. Rix & Zaritsky ( 1995| )), our standard 
mass model has rather narrow and weak spiral arms. This 
might be simply because only the positions of the spiral 
arm tangents are constrained by the COBE/DIRBE L-band 
data, while the spiral arm heads near the bar are not well- 
constrained. It is therefore possible that our standard mass 
model underestimates the strength of the Galactic spiral 
arms. We have therefore investigated a model strongarms 
in which we enlarged non-axisymmetric forces in the disk, 
by multiplying with 1.5 the m > 2 multipoles of the spiral 
arm component (i.e., the component corresponding to the 
density outside of rcut = 3.5 kpc) in the standard model, 
and computed the gas flow in this modified potential. 

We show in Figure ^ the TVC and in Figure ^ the 
(I, v) diagram of this model strongarms. The fit to the TVC 
observations is not as good as for the standard model 20, 
but in the (Z, v) diagram the spiral arm ridges are quite 
similar. Due to the stronger spiral arm gravity, the bumps in 
the terminal velocity are somewhat stronger, but still in the 
acceptable range. Interestingly, the 3kpc-arm now fits the 
observations nearly perfectly. This shows that the strength 
of the spiral arms is an important parameter for the observed 
kinematics of the 3 kpc arm, and the non-circular motions 
in this region of the Galactic (Z, v) diagram cannot simply 
be used to determine the bar aspect angle (cf. Weiner & 
Sellwood 1999). It also suggests that our standard model 
can be improved when a better spiral arm model becomes 
available. 



5.4 Asymmetric Models 

There are indications in the HI and CO surveys that the 
inner Galaxy's gas distribution deviates significantly from 
point-symmetry with respect to the centre. An example is 
the 3kpc-arm, which has no clear counter-arm in the ob- 
served (Z,w) diagrams. That is not to say that there is no 
counter-arm; if asymmetric, its inner parts could for exam- 
ple appear at similar locations in Fig. ^ as the arm which 
reaches the TVC at Z « 25°. All symmetric mass models, 
however, yield counter-arms with about the same absolute 
velocities as the 3kp c-arm . A way out of this dilemma has 
been shown by Fux (1999). In his model, the 3kpc-arm and 



its counter-arm are significantly disturbed by strong non- 
axisymmetric modes in the Galactic centre, as well as in 
the outer disk. Such a mechanism may also help to explain 
that the peak in the terminal velocity at I — +2° appears 
much higher than at I = —2° in the CO data. In addition. 



the X2-disk in the centre may be disturbed by this effect, 
and this might explain the uneven gas distribution seen in 
CS. There is no evidence for an asymmetric mass distribu- 
tion in the NIR data, so in our mass models we can only 
introduce asymmetry in the density by hand and study the 
consequences. 

To this end, we first created uneven m modes in the disk 
by weakening one or two arms in the initial model, or by 
moving two spiral arms closer together. In such gas models 
we observed strong effects on the position of gas shocks and 
their relative strength in the outer disk. Some cases look 
similar to the result of Fux with an almost 3-armed outer 
disk structure. The inner arms, especially the 3kpc-arm, did 
not change much, however. Only in an extreme case have 
we been able to move the 3kpc-arm to higher velocities, 
while simultaneously the counter-arm was moved to lower 
velocities, but this model does not fit the observations well 
overall. 

In a second class of asymmetric models, we let the stel- 
lar bar centre rotate on a circular orbit with radius J?bar 
and pattern speed Qc ~ — 60Gyr~^ = — f2p, i-e., the centre ro- 
tates backwards with respect to the bar and with the same 
pattern speed as the bar. This introduces also a third param- 
eter, the phase Obar of the bar centre rotation at model age 
O.OOGyr. This approach was motivated by A'-body simula- 
tions (Debattista, work in progress). In simulation tumbling- 
bar (fie = -60Gyr"\Rbar = 800pc,abar ~ 80°), we see the 
3kpc-arm and its counter-arm move in the right direction in 
the model's (Z, v) diagram (Figure ^) and the 3 kpc-arm fits 
the observations well. However, we have not found a snap- 
shot of a model at which the 3kpc-arm, its counter-arm and 
the overall spiral pattern all fit the data well. 

Obviously, the available freedom in introducing devia- 
tions from point-symmetry is very large. We have only tested 
a few attractive possibilities, and these models show that 
asymmetries in the Galactic mass distribution may be im- 
portant and could be at the root of some of the remaining 
problems in our standard model. 



6 SUMMARY AND CONCLUSIONS 

We have used new gas flow models to investigate the dynam- 
ics of the Milky Way (MW) Galaxy from observed (Z, v) di- 
agrams. Steady-state gas flows in rotating, point-symmetric 
gravitational potentials for the Galactic bar and disk were 
determined with SPH simulations. The potentials were de- 
rived from non-parametric estimates of the spati al nea r- 
infrared luminosity density (Bissantz & Gerhard (2002)), 



based on the de-re ddened COBE/DIRBE L-band map of 
Spergel et al. (1995), but also in corpo ra ting c lump giant star 
count data from Stanek et al. (1994), (1997). The luminos- 
ity models contain a spiral arm model for the disk, and were 
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Figure 22. (/, v) diagram for the best evolutionary age (0.32Gyr) of model tumblingbar. In this model the centre of the stellar bar is offset 
from the galaxy centre by 300 pc, and rotates around it with —Hp. Note that the 3kpc-arm is reproduced well, but not its counter-arm. 



converted to mass models assuming constant mass-to-light 
ratio in the inner MW. 

Our best gas flow model gives a very good fit to the 
Galactic terminal velocity curve for \l\ > 15°, and to the 
spiral arm ridges in the observed CO (/, v) diagram. This 
has enabled us to investigate a number of dynamically im- 
portant parameters such as the bar and spiral arm pattern 
speeds, the multiplicity of the spiral structure in the poten- 
tial, and the bar angle. The main results from this study are 
as follows. 

1) In gas flow models with separate pattern speeds S7p 
for the bulge/bar and fi^p for the spiral pattern, the spiral 
arms go through the bar corotation region. Thus {I, v) di- 
agrams of such models show well-defined spiral arm shocks 
(ridges) through corotation, next to areas which appear to 
be nearly void of gas. By contrast, in single pattern speed 
models the spiral arms dissolve in the bar corotation region, 
so that the gas fills this region of the {l,v) diagram approx- 
imately evenly, and no voids exist. 

2) Similar voids are visible in the observed ^^CO {l,v) 
diagram. From a comparison with model {l,v) diagrams 
with different flap but similar flp we find evidence for sep- 
arate pattern speeds in the Milky Way. The existence of 
self-consistent models with separate bar and spiral arin_gat- 
tern speeds was demonstrated by Rautiainen & Salo ( 199!: ) 
in a study of two-dimensional N-body simulations, some of 
which included a massless, dissipative gas component. Mod- 
els with a growing bar amplitude also suppor t spir al arms 
in the corotation region (Thielheim & Wolff (1982)); how- 
ever, it is likely that when self-gravity is included and the 
spiral arm amplitude becomes non-linear, the growing spiral 
pattern will again develop an independent pattern speed. 

3) From a series of models the preferred range for the 
bar pattern speed in the MW is Qp = 60 ± 5Gyr"\ cor- 
responding to corotation at 3.4 ± 0.3 kpc. This agrees well 
with previou s pat tern speed deter minations by Englmaier 
fc Ge rhard ( [l999[ ), Dehnen ( ^000| ), and Debattista et al. 
(2002). The bar pattern speed is well constrained because it 
influences not only the inner spiral structure, but also the 
position of two outer spiral arms in the Iv-plot. Models with 
ilp = 50Gyr"^ and ilp = 70Gyr"^ are inferior. 

The spiral arm pattern speed is less well constrained. 
Our preferred value is Qap ~ 20Gyr"^, but models with 
larger Q^p < ^p give only marginally inferior flts to the ob- 
served {I, v) diagram. 

4) Gas flows in models which include massive spiral 
arms clearly fit the observed ^^CO {l,v) plot better than if 
the potential does not include spiral structure. Furthermore, 
comparing models with two and four arms in the gravita- 
tional potential, we found that only models with four mas- 
sive arms reproduce the Galactic {l,v) diagram, while gas 
flows in two-armed potentials do not resemble the spiral arm 
pattern of the Milky Way. 

In Galactic models with four-armed potentials and sep- 
arate spiral arm pattern speed, the gas flow has two pairs 
of inner arms which rotate with the bar (lateral, and corre- 
sponding to the 3 kpc arm), and four outer spiral arms which 
exhibit a complicated, time-dependent back-and-forth oscil- 



lation in the bar frame. The outer and inner spiral structures 
are connected by a time-dependent transition region around 
bar corotation. 

5) From a further series of gas models computed for 
different bar angles, using separately determined luminos- 
ity mode ls and gravitational potentials as in Bissantz & 
Gerhard (2002), we found a range of acceptable bar angles 
20° < (^bar < 25°. The models for (pta, = 15° and 30° are 
clearly inferior, which is mainly due to differences in the 
inferred gravitational potential. 

The model with ((^ta, = 20°, D.p = GOGyr^S D,sp = 
20Gyr~^) gives an excellent flt to the Galactic terminal ve- 
locity curve for 10° < |/| < 50°, and to the gaps and spiral 
arm ridges in the observed CO (/, v) diagram. There are still 
discrepancies in the bar corotation region where the poten- 
tial has uncertainties: the 3 kpc arm has too low non-circular 
velocity, and its counterarm is missing in the data. In the 
bulge region, closed orbits reproduce the TVC well, while 
the gas model has lower velocities. This may be a resolution 
problem in the SPH model, but could in part also be due 
to uncertainties in the potential which influence the orbit 
shapes. 

6) The 3 kpc arm non-circular velocities can be re- 
produced by a model in which we artiflcally increased the 
m > 2-multipoles of the spiral potential component by 1.5 
while keeping all other dynamical parameters fixed. This is 
well within the uncertainties. Guided by a number of asym- 
metries in the observed Milky Way gas distribution, we also 
investigated potentials which are no longer point-symmetric. 
Some of these models improved the fit to the 3kpc-arm and 
its counter-arm. Although we did not find a model which at 
the same time reproduces the entire {I, v) diagram as well as 
our standard model, these models suggest that such asym- 
metries may be important for better understanding the gas 
flow in the inner Milky Way. 
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